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The  pressure  induced  infiltration  of  nanoporous  materials  has  recently  received  increasing 
attention.  When  hydrophobic  nanoporous  materials  are  immersed  in  water,  at  the  atmosphere  pressure, 
pat,  due  to  the  capillary  effect  soaking  does  not  occur  spontaneously.  As  the  pressure  increases  to  a 
critical  value,  pm,  water  can  infiltrate  into  the  nanopores,  accompanied  by  a  large  increase  in  solid-liquid 
interfacial  energy.  When  the  pressure  is  reduced  back  to  pat,  however,  the  water  molecules  may  not  come 
out  of  the  nanoenvironment,  depending  on  the  factors  that  are  still  under  investigation.  According  to  the 
limited  data  currently  available  in  open  literature,  the  “nonoutflow"  is  likely  to  occur  when  the  pore  size  is 
in  the  mesoporous  range  (2-80  nm),  and  the  “outflow”  is  relatively  easy  if  the  pore  size  is  smaller  than 
about  1-2  nm. 

The  conventional  interface  and  microfluidic  theories,  such  as  the  Washburn  type  analyses  in 
which  the  solid-liquid  interfaces  are  assumed  to  be  isolated,  have  failed  in  explaining  this  size  effect,  as 
well  as  the  thermal  effect,  the  recovery  behavior,  the  dynamic  behavior,  the  system  selectivity,  and  the 
gas  diffusion  phenomenon  that  will  be  discussed  shortly.  In  order  to  model  the  pressure  induced 
infiltration,  a  few  frameworks,  such  as  the  evolution  of  pore  clusters,  the  flow-direction  dependence  of 
contact  angle,  and  the  phase  transformation  with  constant  system  volume  and  zero  gas-phase  nucleation 
barrier,  have  been  proposed.  However,  currently,  since  the  experimental  setups  are  usually  sophisticated 
and  the  direct  observation  of  infiltration  processes  is  difficult,  there  are  few  experimental  evidences  that 
can  support  these  theories.  For  instance,  at  the  nm  scale,  the  basis  of  continuum  fluid  mechanics 
involved  in  the  effective  phase  transformation  and  the  flow  direction  theory  may  no  longer  be  valid,  and 
ignoring  the  system  volume  change  and  the  energy  exchange  between  gas/liquid  phases  can  be 
questionable. 

We  designed  a  simple  transparent  poly(methyl  methacrylate)  (PMMA)  system  so  as  to  directly 
observe  the  infiltration/defiltration  behaviors  of  nanoporous  particles.  The  material  under  investigation 
was  end-capped  Fluka  100  C8  reversed  phase  mesoporous  silica  with  the  average  pore  size  7  =  7.8  nm 
and  the  standard  deviation  6 r  =  2.4  nm.  The  specific  pore  volume  was  560  mm3/g,  and  the  specific 
surface  area  was  287  m2/g.  The  Barrett-Joyner-Halenda  (BJH)  adsorption  characterization  measurement 
was  performed  at  the  Quantachrome  Instruments.  The  surface  coverage  was  10-12%  (±  4  pmol/m2), 
which  led  to  a  high  degree  of  hydrophobicity.  The  particle  size  was  in  the  range  of  15-35  pm.  Prior  to  the 
infiltration  tests,  the  silica  particles  had  been  heated  in  air  at  150°C  for  12  hours. 

The  aqueous  suspension  of  0.5  g  of  the  mesoporous  silica  particles  was  sealed  in  the  PMMA 
cylinder  by  a  stainless  steel  piston  with  reinforced  gaskets.  Initially,  no  air  bubble  could  be  observed.  The 
infiltration  experiment  was  performed  using  an  Instron  5569  machine.  The  piston  was  first  compressed 
into  the  container  at  a  constant  rate  of  1.0  mm/min.  Once  the  pressure  exceeded  about  50  MPa,  the 
crosshead  was  moved  back  at  the  same  speed.  The  loading-unloading  cycle  was  repeated  until  the 
absorption  isotherm  curves  converged  to  the  steady-state,  as  shown  by  curves  (a,  b)  in  Fig.  1 ,  where  the 
specific  volume  variation  is  defined  as  A  Vo IW,  with  A  V0  being  the  volume  change  of  the  system  and  Wthe 
mass  of  the  silica  gel.  The  system  was  then  thermally  treated  in  a  temperature  bath  in  the  range  of  30- 
80°C  for  0.5  hour  either  (1)  immediately  or  (2)  after  resting  at  room  temperature  for  6-24  hours,  followed 
by  another  loading-unloading  test.  The  testing  results  are  shown  by  curves  (c,  d)  in  Fig.  1 ,  respectively. 
Altogether  four  samples  were  tested. 

Curve  (a)  in  Fig.  1  shows  that,  following  the  initial  linear  stage,  as  the  pressure  reached  the 
infiltration  pressure,  p,„  =  17  MPa,  the  water  was  forced  into  the  nanopores,  causing  the  large  increase  in 
system  compressibility.  If  the  pore  size  were  perfectly  uniform,  the  plateau  should  be  flat.  In  the  current 
system,  due  to  the  pore  size  distribution,  the  slope  of  the  absorption  isotherm  was  finite.  Eventually,  at 
about  30  MPa,  most  of  the  pores  were  filled  and  the  system  compressibility  decreased  rapidly.  The 
volume  variation  associated  with  the  plateau  region  was  about  0.55  cm3/g,  close  to  the  BJH  result  of  the 
specific  pore  volume.  As  the  pressure  was  reduced,  the  confined  water  remained  in  the  nanopores,  and 
thus  the  unloading  curves  were  quite  linear.  During  the  loading-unloading  process,  there  was  no 
significant  change  in  system  appearance,  except  for  the  variation  in  volume,  indicating  that  the  gas 
entrapped  in  the  nanopores  dissolved  in  the  liquid  phase.  Since  no  air  bubble  could  be  observed  even 
when  the  pressure  was  reduced  back  to  the  atmosphere  pressure,  the  gas  content  in  the  liquid  phase 
outside  the  nanoporous  particles  must  be  quite  constant;  that  is,  the  gas  molecules  remained  in  the 
nanopores.  Because  of  the  "nonoutflow”,  the  extent  of  infiltration  was  considerably  lowered  in  the 
following  loading-unloading  cycles  (see  curve  b  in  Fig.1). 
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When,  immediately  after  the  first  loading-unloading  cycle,  the  system  was  thermally  treated,  the 
liquid  phase  was  still  clear  and  little  air  bubbles  were  formed.  However,  the  energy  absorption  capacity  of 
the  system  was  recovered  significantly  (see  curve  c  in  Fig.  1  and  Table  I),  indicating  clearly  that  after  the 
thermal  treatment  a  certain  portion  of  the  porous  space  was  occupied  by  gas  phase.  When  the 
temperature  exceeded  50°C,  the  system  could  be  almost  fully  recovered. 

If,  on  the  other  hand,  after  the  first  infiltration  cycle,  the  specimen  was  rested  at  room  temperature 
under  pat,  the  system  appearance  would  change  gradually.  After  6  hours,  a  large  number  of  air  bubbles 
with  the  sizes  in  the  range  of  0.05  mm  to  0.5  mm  were  formed  and  therefore  the  sample  was  no  longer 
transparent.  After  24  hours,  the  total  volume  of  the  air  bubbles  was  estimated  as  0.25  cm3/g,  while  the 
energy  absorption  capacity  was  still  close  to  zero,  suggesting  that  the  nanopores  were  still  filled  by  liquid. 
Thermally  treating  this  system  would  cause  a  partial  recovery,  as  shown  by  curve  (d)  in  Fig.  1 ,  where  the 
system  recoverability,  Rp,  is  defined  as  EJEU  with  Et  being  the  absorbed  energy  in  the  loading¬ 
unloading  cycle  subsequent  to  the  postponed  thermal  treatment;  the  treatment  temperature  was  set  to 
70°C;  and  t,  is  the  room-temperature  resting  time.  Clearly,  during  the  room-temperature  resting,  a  certain 
amount  of  gas  molecules  diffused  out  of  the  nanopores,  and  the  decrease  in  system  recoverability  should 
be  attributed  to  the  reduced  excess  gas  content  in  the  nanoenvironment. 

The  motion  of  the  gas-liquid  contact  line  in  a  nanopore  can  be  considered  as  the  result  of 
expansion  or  shrinkage  of  the  gas  phase.  Note  that  at  the  nm  level,  sharp  liquid-gas  interfaces  may  not 
exist.  Nevertheless,  effective  boundaries  can  be  defined,  e.g.  in  the  context  of  Gibbs  dividing  surfaces,  as 
shown  in  Fig.2.  In  an  initially  filled  nanopore,  the  formation  of  a  gas  phase  nucleus  increases  the  free 
energy  of  system  by  (a)  A p-  VG,  where  A p  is  the  specific  nucleation  energy  and  VQ  =  n?h  is  the  volume  of 
the  gas  phase,  with  h  being  the  nucleus  length;  (b)  jqL  AGl,  where  7feL  is  the  surface  tension  of  the  liquid 
and  Agl  =  is  the  gas-liquid  interface  area;  and  (c)  the  external  work  p AV,  where  p  is  the  applied 
pressure  and  AV  =  VG  is  the  system  volume  change.  On  the  other  hand,  with  the  nucleation  of  the  gas 
phase,  since  the  liquid  is  nonwetting,  the  liquid-solid  interfacial  energy  is  reduced  by  A/Aqs,  where  A/is 
the  difference  between  the  gas-solid  and  liquid-solid  interfacial  energy,  and  AGs  =  2 nrh  is  the  gas-solid 
interface  area.  Note  that,  if  the  effective  gas-liquid  interfaces  are  irregular,  geometry  factors  should  be 
used  in  the  calculations  of  AGl  and  VG,  which  are  ignored  in  the  following  discussion  for  the  sake  of 
simplicity. 

The  thermodynamics  equilibrium  condition  can  then  be  stated  as 

ApVG+fciAGim  +  pAV=A},AGs,  (1) 

which  can  be  rewritten  as 

r  =  ra  (2) 

or 

h  =  ha  (3) 

where  r  - _ ^Ay  _ is  the  critical  pore  radius  and  /,  -  ^cl  is  the  critical  nucleus 

(p  +  Ap)  +  2yGL/h  "  lAy/r-ip  +  Ap) 

size  of  gas  phase.  According  to  Eq.(2),  when  r  <  rc„  the  gas  phase  is  stable  and  can  expand 
spontaneously,  which  eventually  leads  to  the  “outflow”;  when  r  >  rct,  the  liquid-to-gas  phase 
transformation  is  energetically  unfavorable  and  therefore  the  gas  phase  nucleus  will  vanish,  resulting  in 
the  “nonoutflow”.  The  relationship  between  rci  and  h  is  shown  in  Fig.3.  The  value  of  A?' is  estimated  as 
pmr  / 2 ,  and  Ap  is  taken  as  po  pq,  where  po  =  2257.1  J/kg  is  the  specific  energy  of  evaporation  of  water 

and  pa  is  the  density  of  gas  phase,  which,  if  we  assume  that  the  gas  phase  follows  the  law  of  ideal  gas, 
can  be  assessed  as  pc  =  pb/(1+p/Pat).  with  po  =  1.21  kg/m3  being  the  air  density  at  pat.  For  the  first-order 
approximation,  the  macroscopic  value  of  fcL,  72  mJ/m2,  is  used.  Figure  3  indicates  that  decreasing  p  or 
increasing  h  is  beneficial  to  keeping  r  <  rc,  .  According  to  Eq.(3),  there  exists  an  ultimate  pore  size 

r  (  -  — X.  ,  at  which  hCI  — »  » .  if  r  >  ruh  the  “outflow”  is  impossible  for  any  [Ay  Ap,  jfeL,  p}. 
p  +  Ap 

Based  on  a  Fourier  Transform  Infrared  Spectroscopy  (FTIR)  analysis,  it  was  confirmed  that,  at 
liquid-solid  interfaces,  there  exist  a  large  number  of  nm-scale  gas  phase  nuclei.11  As  temperature  rises, 
the  fraction  of  bigger  nuclei  increases,  i.e.  the  average  h  is  larger,  which  explains  that,  under  pat,  while  at 
room  temperature  the  “outflow”  was  negligible,  when  T >  50°C  the  gas  phase  nucleation  and  growth  could 
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occur.  After  the  room-temperature  resting,  however,  due  to  the  significant  decrease  in  gas  content  in  the 
nanopores,  the  average  size  of  gas  phase  nuclei  is  smaller,  and  therefore  the  same  thermal  treatment 
can  cause  only  a  reduced  system  recovery. 

Future  Work 

The  proposed  research  and  work  plan  are  not  affected  by  the  transfer  of  the  place  of  operation. 
The  objective  of  the  proposed  research  program  is  to  provide  a  scientific  basis  for  the  development  of 
high-performance  nanoporous  volume-memory  systems  that  are  of  superior  energy  density  and 
deformability.  The  effects  of  admixtures,  surface  properties,  porous  structures,  temperature,  and  potential 
difference  on  the  solid-liquid  interfacial  energy  and  the  kinetics  of  infiltration  will  be  investigated 
systematically  through  experiments  and  computer  simulations. 

In  the  first  period,  we  will  focus  on: 

(1)  Understanding  the  behaviors  of  nanoporous  carbon  based  systems  through  infiltration 
measurements.  Nanoporous  carbon  particles  and  membranes  with  various  pore  sizes  and  pore 
volume  fractions  will  be  tested. 

(2)  Performing  infiltration  experiment  on  nanoporous  zeolite  and  silica  based  systems  at  various 
temperatures,  through  which  the  thermal  effect  can  be  analyzed. 

(3)  Performing  staged  infiltration  experiment  on  nanoporous  systems  to  analyze  the  kinetics  of  infiltration 
under  various  conditions. 

In  the  second  period,  based  on  the  experimental  results  obtained  in  year  1 ,  we  will  further  investigate  the 
nanoporous  intelligent  systems: 

(1)  Continue  the  infiltration  experiment  on  nanoporous  systems  to  evaluate  different  surface  treatment 
techniques  and  admixtures. 

(2)  Continue  the  staged  infiltration  experiment  to  analyze  the  “flow”  rate  in  nanopores,  with  the  pore  size 
and  surface  species  varying  in  broad  ranges.  The  effects  of  the  properties  of  gas  species  will  be 
investigated. 

(3)  Evaluate  the  effectiveness  of  intrinsic  and  extrinsic  factors  through  atomic  simulations.  The  results  will 
be  compared  with  the  experimental  data. 

(4)  Perform  numerical  analyses  on  the  evolution  of  effective  pore  clusters,  such  that  the  atomic  behaviors 
in  nanopores  can  be  related  to  the  aggregate  response  of  nanoporous  particles. 

In  the  third  period,  with  the  important  factors  being  fundamentally  much  better  understood,  the  study  will 
lead  to  the  eventual  optimization  of  the  systems: 

(1)  Complete  the  experiment  on  intelligent  systems  of  different  admixtures  and  surface  properties. 

(2)  Develop  the  multiscale  model. 

(3)  Complete  the  parameterized  study  and  identify  the  optimum  materials  and  processing  conditions. 

(4)  If  necessary,  investigate  nanoporous  systems  based  on  metal  doped  zeolites  so  as  to  control  the 
conductivity  of  electrodes. 

(5)  Complete  reporting  of  experimental  and  simulation  results. 
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Table  I.  The  system  recoverability,  Rs=  E 2IEi,  as  a  function  of  the  thermal  treatment  temperature,  T, 
where  E  ,•  is  the  absorbed  energy  in  the  /'- th  loading-unloading  cycle  (/'  =  1 ,2). 
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Fig.1  The  sorption  isotherm  curves:  (a)  the  first  loading-unloading  cycle  (the  dashed  line);  (b) 
the  second  (the  dotted  line),  the  third,  and  the  fourth  loading-unloading  cycles  without 
thermal  treatment;  (c)  after  the  immediate  thermal  treatment;  and  (d)  after  the  postponed 
thermal  treatment. 
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Fig.2  A  schematic  diagram  of  the  confined  phase  transformation  in  a  nanopore. 
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